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Abstract: Finite ’t Hooft coupling corrections to multiple physical observables in 
strongly coupled A/” = 4 supersymmetric Yang-Mills plasma are examined, in an attempt 
to assess the stability of the expansion in inverse powers of the ’t Hooft coupling A. 
Observables considered include thermodynamic quantities, transport coefficients, and 
quasinormal mode frequencies. Although large A expansions for quasinormal mode 
frequencies are notably less well behaved than the expansions of other quantities, we 
find that a partial resummation of higher order corrections can significantly reduce the 
sensitivity of the results to the value of A. 
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1 Introduction 

Relativistic heavy ion collisions at the LHC and RHIC colliders produce a novel state 
of matter, the quark gluon plasma (QGP) [1, 2]. At accessible energies, the produced 
plasma is strongly coupled — so strongly coupled that it cannot be described in terms 
of long-lived quasiparticles. Evidence for this comes from the effectiveness of hydrody¬ 
namics in modeling experimental results (for a recent review, see ref. [3]) with viscosity 
close to the holographic value [4], as well as from the very short values of screening 
lengths in hot QCD which are computable using lattice gauge theory [5]. 

Hydrodynamics provides an effective description of the dynamics of the plasma at 
sufficiently late times after the collision, but is not applicable to early time dynamics 
when the produced plasma is very far from local equilibrium. Moreover, hydrodynamics 
is not adequate for understanding the important physics of hard probes of the medium. 
Unfortunately, alternative theoretical approaches for calculating, reliably, properties of 
a strongly coupled plasma are very limited. 

In recent years, gauge/gravity duality (or “holography”) has provided a new tool 
for understanding strongly coupled systems. In its simplest and most studied form, 
gauge/gravity duality relates properties of maximally supersymmetric SU{Nc) Yang- 
Mills theory (A/* = 4 SYM), in the A), —)■ oo limit, to gravitational dynamics of higher 
dimensional asymptotically anti-de Sitter spacetimes. Under this duality, the process 
of equilibration and thermalization in the quantum held theory is precisely related to 
gravitational dynamics involving the formation and subsequent equilibration of black 
hole horizons. 
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Much work has been done using gauge/gravity duality to study aspects of strongly 
coupled dynamics relevant to heavy ion collisions; see, for example, refs. [6-8] and 
references therein. This includes calculations of the drag on a heavy quark [9-13] or 
light quark [14, 15] propagating through a strongly coupled medium, jet quenching 
[16, 17], particle production [18-22], isotropization dynamics [23-26], boost-invariant 
flow [27, 28], collapsing bulk scalar helds, planar shells, and balls of dust [29-33], 
collisions of planar shock waves [34-37], and collisions of fully localized shock waves 
resembling Lorentz contracted nuclei [38, 39]. Because of the precise mapping between 
gauge and gravitational dynamics provided by gauge/gravity duality, in all this work 
one is honestly computing properties of a strongly coupled non-Abelian gauge theory. 
There is just one problem — the theory in which these calculations are done, namely 
A7 = 4 SYM, is not real QCD. 

The strong coupling, large Nc limit of A/' = 4 SYM, to which gauge/gravity duality 
applies, may be viewed as a three-step deformation of QCD: (i) the fundamental repre¬ 
sentation quarks of QCD are replaced by a collection of adjoint representation matter 
helds, both fermions and scalars, thereby turning QCD into A/" = 4 SYM, (ii) the ’t 
Hooft coupling A = which no longer runs with energy scale in TV = 4 SYM, is 

tuned to very large values, and (hi) the gauge group rank, Nc, is sent to inhnity. Qual¬ 
itative properties of the deconhned plasma phase are stable under these deformations: 
the high temperature phase of the theory remains a non-Abelian plasma with Debye 
screening and a hnite correlation length; spacelike Wilson loops continue to show area 
law behavior; and long distance, low frequency dynamics continues to be described by 
neutral huid hydrodynamics. 

Lattice gauge theory simulations have shown that thermodynamic properties of 
SU{Nc) Yang-Mills plasma scale very smoothly with Nc [40, 41], suggesting that the 
large Nc limit should be well-behaved for most observables of interest, and moreover 
that the SU{3) theory is already fairly close to the Nc = oo asymptotic limit. Where 
results from hot QCD lattice simulations (in the experimentally relevant temperature 
range, 1.5 < T/Tc < 4) are available to be compared to holographic computations in 
A7 = 4 SYM, a variety of important physical quantities such as the equation of state, 
ratios of screening masses to temperature in various symmetry channels, and estimates 
of the shear viscosity to entropy density ratio, rj/s, show agreement to within at least 
a factor of two, and often much better. 

Consequently, in the above deformations which connect holographic models to 
QCD, the step which likely produces the largest changes in thermal properties, and 
about which the least is known, is step (ii): sending the ’t Hooft coupling to values 
large compared to unity. At the relevant energy scales in hot QCD, the appropriate 
value of the ’t Hooft coupling (in physically sensible schemes) is presumably some- 
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where in the range 10-40 — corresponding to Og = between 0.3 and 1, not 

some truly enormous number. Therefore, improved understanding of the dependence 
of physical quantities in A/" = 4 SYM on the value of A is highly desirable. It is known 
that hnite A corrections appear in the form of inverse fractional powers, beginning with 
In this paper, we collect, extend, and examine available results for £nite-A cor¬ 
rections to thermal observables in an effort to gain some insight into the stability of 
the expansion in inverse powers of A and the applicability of holographic predictions to 
physics at realistic values of the’t Hooft coupling. 

The paper is organized as follows: In section 2, we summarize and discuss hrst order 
£nite-A corrections to a variety of thermal observables. A basic observation is that the 
relative size of the hrst £nite-A correction is substantially larger for quasinormal mode 
(QNM) frequencies than for other observables. Section 3 then recaps the holographic 
calculation of two point correlation functions, from which transport coefficients and 
quasinormal mode frequencies are extracted, in a manner which allows one to extract 
the hrst order hnite-A correction or perform a partial resummation of higher order 
hnite-A corrections. Results of the two procedures are shown, and compared, for the 
hrst few quasinormal modes of the current-current correlator and the shear channel of 
the stress-energy correlator, as well as for the plasma conductivity and shear viscosity. 
The hnal section 4 contains a few concluding remarks. 

2 First order corrections: collected results 

Considerable prior work exists examining hnite-A corrections to holographic results. 
This includes analyses of the equation of state [42], shear viscosity 7] [43, 44], plasma 
conductivity a [45], photon production and transport [46, 47], and various higher-order 
transport coefficients [48]. More recent work has also considered hnite-A corrections to 
quasinormal mode frequencies and oh-equilibrium spectral densities obtained from the 
current-current and stress-energy correlators [49-51]. 

The hnite coupling corrections appear as a power series in with the hrst 

corrections being proportional to It proves convenient to dehne the constant 

7 = |C(3)-^ = 8 ^(3) (s'ym-^c) , (2.1) 

where |C(3) ~ 0.15. The benchmark range of 10-40 for the ’t Hooft coupling A 
corresponds to values of 7 between about 0.005 and 0.0006. 

In table 1, we collect the values of the hrst two terms in the strong coupling 
expansions of various thermal observables: the entropy density s, the shear viscosity 
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Quantity 

0(7°) 

0(7Q 

Reference 


1 

15 7 

|42] 


1 

135 7 

[44] 

An rj/s 

1 

120 7 

[44] 

ailaEuN^T)-^ 

1 

14993/9 7 

[45] 

Fo (ctEMiV'T^)-' 

0.053678 

23.53797 

This work 

Fi (aEuN^T^)-^ 

0.472771 

-224.4698 7 

This work 

^hear (27rT)-^ 

2.585 - 2.382 i 

(1.029 +0.9571) 10^7 

|50] 

(U™ (27rT)-^ 

2-2i 

(1.34+ 0.43Q 10^7 

|51] 


Table 1. Zeroth and first order terms in the expansion of various thermal observables in 
powers of 7 = Results are shown for the entropy density s, shear viscosity 

7 , viscosity to entropy density ratio rj/s, electrical conductivity a, the first two moments, 
To and Ti, of the photoemission spectrum, and the second quasinormal mode frequencies, 
cj™ and at zero wavevector, for the electromagnetic current and shear channel of the 

stress-energy correlator, respectively. 

r] and viscosity to entropy density ratio rj/s, the electrical conductivity a,^ and the 
second quasinormal mode frequencies ca™ and for the electromagnetic current 

and the shear channel of the stress-energy correlator, respectively. Also included are 
results for the hrst two moments of the photoemission spectrum, dehned as 

r°° r/r 

( 2 . 2 ) 

Jo dk 

For the entropy density (or equivalently, the pressure p = sT/4), the hrst hnite- 
coupling correction is modest; the 0(7^) term does not exceed the leading 0 ( 7 *^) term 
as long as A > 1.72. For the shear viscosity or viscosity ratio rj/s, the corresponding 

^To define the SYM electromagnetic current and associated electrical conductivity, one weakly 
gaugesa?7(l) subgroup of the global S'[/(4)fl; flavor symmetry ofA/’ = 4SYM [18]; see the next section 
for details. 
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crossover points where the first corrections equal the leading term occur at A ~ 7.4 or 
6.9, respectively. For the electric conductivity, this crossover lies at A ~ 39.7, while the 
crossovers for the photoemission moments Fq and Fi are at 16.3 and 17.2, respectively. 
All these values are below, or at least within, our 10-40 range of benchmark values for A. 
The situation, however, is rather different for the quasinormal mode frequencies shown 
in table 1. The hrst order corrections exceed the leading order term when A < 71.2 
(for or A < 382.2 (for suggesting that their A = cx) limits are likely to give 

poor predictions for the values of these quantities in the phenomenologically interesting 
range of’t Hooft couplings. 

A priori, it is not clear whether the above behavior of the QNM frequencies is due 
to an abnormally large first term in an otherwise well-behaved expansion, or whether 
the quantities in question are particularly sensitive to hnite coupling corrections, so 
that their expansions in 7 have an abnormally small range of utility. Deciding between 
these alternatives would, in principle, require an all orders determination of the strong 
coupling expansion, which is far beyond the reach of present day technology. In this 
paper, we adopt a far more modest goal. We will investigate a simple resummation 
applicable to quasinormal mode frequencies and related observables which takes into 
account a subset of higher order terms in the expansion in powers of 7 , and see if 
this improves the behavior of the resulting series. At the very least, this investigation 
should be helpful in inferring the range of utility of the above first order results. 

3 Finite coupling corrections to correlators and QNMs 

Finite coupling corrections to thermal observables are generated by higher derivative 
(or a') corrections to the 10-dimensional type IIB supergravity action, which takes the 
schematic form 

SiiB = S<IB + 7 SiiB + 7*''= dm ' + • • • ■ (3-1) 

The first order correction S'jjg includes fourth powers of the Riemann tensor plus 
terms, related by supersymmetry, that involve the self-dual five form (see, for example, 
ref. [52]). For the free energy (or pressure), it is sufficient to evaluate the first order 
correction terms in the action on the unmodihed AdS-Schwarzschild solution [42]. For 
other observables, one must insert into the lOD action (3.1) an appropriate ansatz for 
the lOD metric, hve form held strength, and any other helds relevant for the observ¬ 
able of interest. A Kaluza-Klein reduction eliminating the compact internal space (for 
physics which only depends on the lowest KK modes) leads to a 5D a'-corrected action 
for the relevant bulk helds in asymptotically AdS spacetime. Using the corrected 5D 
action, observables of interest are computed using the standard holographic correspon¬ 
dence. This typically involves deriving a' corrected equations of motion for the relevant 
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bulk fields and then solving these equations order by order in 7 . For more details of 
such calculations see, for example, refs. [42, 44, 45, 51, 53]. 

Below, we illustrate explicitly the above procedure as applied to the calculation of 
the electromagnetic current and stress-energy correlators and the extraction of their 
associated QNM spectra and transport coefficients. In each case, we first perform the 
computation in a way that consistently truncates the result after the linear term in 7 . 
Thereafter, we present an alternative calculational scheme which resums a subset of 
higher order corrections, all originating from the correction term to the supergrav¬ 
ity action. We emphasize that, as the explicit forms (and physical effects) of higher 
order terms in the supergravity action are presently unknown, our resummation only 
captures a limited subset of corrections involving higher powers of 7 . Nevertheless, 
the results of this partial resummation will be seen to have interesting and suggestive 
implications for the stability of holographic predictions at phenomenologically relevant 
values of the ‘t Hooft coupling. 


3.1 Current-current correlator 


We begin by considering correlators of the electromagnetic current operator of 
AT = 4 SYM, defined by gauging a U{1) subgroup of the S'17(4 )r flavor symmetry.^ 
This current is dual to a U{1) vector field Am in the gravitational description. To com¬ 
pute the two-point correlator (jEMlEM)) must solve for the behavior of linearized 
fluctuations of this bulk gauge field in the background geometry corresponding to the 
equilibrium state of interest. In the near boundary expansion of the bulk gauge field 
Am, the coefficient of the leading term represents a source coupled to the conserved cur¬ 
rent jgjyj, and the coefficient of the first sub leading term encodes the expectation value 
of the current in the presence of this source. Hence, the two-point correlator is given 
by the variation of the subleading coefficient with respect to the leading coefficient. 

In the A —>■ cx) limit, the action for the bulk gauge field is just the standard Maxwell 
action, which leads to the usual (curved space) Maxwell equation. 




(3.2) 


The Af = 4 thermal equilibrium state (in the absence of chemical potentials) is dual to 
the AdS-Schwarzschild geometry, whose metric may be written in the form 

2 T 2 

ds^ = \-f{u) dt^ dy'^ dz'^] + du^. (3.3) 

L^u ^ 4 m ^/( m ) 

^Specifically, we choose the 17(1) subgroup for which the ^7 = 4 SYM fermions have charges 
{!,—1,0, 0}; the explicit form of the current in terms of the ^7 = 4 SYM fields is shown in eq. (2.1) 
of ref. [18]. 
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Here, L is the AdS curvature scale, which we choose to set to unity, while the non- 
extremality parameter is related to the held theory temperature T hy rh = nTL'^. The 
coordinate u is an inverted radial coordinate, with the spacetime boundary lying at 
M = 0; in terms of a conventional non-inverted radial coordinate r, n = The 

blackening function hnally reads /(n) = 1—v?, and vanishes at the black brane horizon 
located at n = 1. 

When evaluated in the above geometry and Fourier transformed with respect to 
the boundary (Minkowski) coordinates, the Maxwell equation (3.2) reduces to a pair 
of decoupled linear ordinary differential equations for the longitudinal and transverse 
components of the electric held [18], 


0 

0 


A_l(m) + A_l(m) H-— A 


A[|'(n) + 


fiu) 

/(n)(w2 -g2) 


uf{uY 

- q^f{u) 


E[Au) + 


uf{u) 


Ell 


(3.4) 

(3.5) 


with (h = oj/{2'kT) and q = q/{27iT) denoting the rescaled frequency and spatial 
wavevector, respectively. Focusing on the transverse electric held (which determines the 
transverse part of the current-current correlator and thus the photoemission spectrum), 
the simple held redehnition \k(M) = f(u) E±(u) converts eq. (3.4) into a Schrodinger- 
like equation at zero energy. 


- T"(n) + l/(n) T(n) = 0, 


(3.6) 


with 


V(u) 


u + Cj"^ — (ff{u) 
uf{uY 


(3.7) 


At non-zero temperature, the retarded current-current correlator may be decom¬ 
posed into transverse and longitudinal pieces via 


=7^i(w,q)n±((h,g) + Pj^(a;,q)n||(a;,g), (3.8) 

with the symmetric projectors dehned by P(^(a;, q) = 0, P^ioj, q) = 5ij — qiqj/q^, and 
q) = - Q^QuIQ"^ - q), where Q = (w, q) and -F q^. The 

transverse correlation function is then given by [54] 


n±(d/,g) 


lim 


E\ 


u) 


u^o E± (u) 


lim 


T'(m) 


u^O \ 1 /( m ) 


(3.9) 


where E± (or T) is the solution to eq. (3.4) [or (3.6)] satisfying infalling boundary 
conditions at the horizon. 
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Numerous physical observables of interest can be extracted from the current-current 
correlator. Quasinormal modes are poles of the retarded correlator, regarded as func¬ 
tions of the complex frequency u for hxed wavevector q, while pole positions at imag¬ 
inary values of q and a; = 0 give thermal screening masses [55]. The zero-frequency 
slope, at vanishing wavevector, determines the electric conductivity. 


a = — lim Im — g=0) 

a;—^-O OJ 


e^T 

Ibvr 


lim lim — 

UJ —^0 —^0 (jlJ 


Im 


T'(n) 


(3.10) 


where e is the (arbitrarily weak) coupling constant of the electromagnetic U{1) gauge 
held coupled to the conserved current.^ Finally, the (equilibrium) photoemission spec¬ 
trum is determined by the imaginary part of the transverse correlator evaluated on the 
lightcone [18], 

^ = ^kn^{k) (-41mnx)|_^^/(2.T)> (3-11) 

where nb(ci;) = — 1)“^ is the usual Bose distribution function. 

Expression (3.9) shows that the correlator will have poles at values of q and u for 
which the denominator, equal to the boundary value of the electric held (or T), vanishes. 
In other words, QNMs represent homogeneous solutions of the bulk Maxwell equations 
satisfying infalling boundary conditions at the horizon and a Dirichlet condition at the 
spacetime boundary. At g = 0, one may solve the transverse equation (3.4) and hnd 
the resulting roots of E± analytically [56] . The result is the famous linear spectrum. 


uJn = n (±1 - i ), 


(3.12) 


while numerical results for non-zero wavevector may be found in ref. [57]. 

To incorporate hnite coupling corrections in the above calculation, one begins with 
the expansion (3.1) of the lOD supergravity action and retains both the leading and 
hrst subleading terms. Schematically, 

dm = ^ V=G [«,„ - . (3.13) 

Sul = 1^ [ (C + T)" , (3.14) 

2^10 J 

where /tio is the lOD Newton constant, Rio the lOD Ricci-scalar, 0 the dilaton held, 
and Fs the hve-form held strength, while C stands for the Weyl tensor, and T for a 
tensor built from the gradient of plus terms quadratic in [42, 52, 53]. After a 

one regards the i7(l) current as a global symmetry current, not coupled to a dynamical electro¬ 
magnetic gauge field, then the associated charge diffusion constant is related to the conductivity by 
the Einstein relation, D = tT/(e^S), where S = | is the A/’ = 4 SYM charge susceptibility. 











reduction to 5D and the extraction of terms at most quadratic in the emergent 5D 
bulk gauge held dual to the EM conserved current, one eventually hnds a 7 -corrected 
Maxwell equation which, once again, may be put into the Schrodinger-like form (3.6) 
[46, 51]. The needed held redehnition becomes T(m) = J2{u)E±{u ), with 


S(f) 


\//('“) 

1 + 7p(n) ’ 


p{u) 


^ [11700 - 343897 - 37760 + 87539 n^] , 

288 

(3.15) 


while the resulting 7 -corrected ehective potential reads 


V(u) = 


u + u 


q^f{u) 


7 


144/(n) 


11700 + 2098482^2 _ 4752055 + 1838319 m® 


+ f (4770 + 11700 m - 953781 m^ + 1011173 m®) 

(4770 M"^ + 28170 m- 1199223 m®) . 


(3.16) 


Given the above potential, solutions to the equation T" = G 4/ may be expanded in 
a power series in 7 . The black brane horizon at m = 1 is a regular singular point of the 
equation, where infalling solutions behave locally as 4/(m) ~ (1— m)*” with characteristic 
exponent r = |(1 — iu). Hence, one may expand the solutions of interest as 

T(m) = (1-m)" [<h(°)(M) + 7<h(®)(M) + •••], (3.17) 

where and have the near-horizon Frobenius expansions, 

00 00 

4>(®)(m) ~ (1-m)", $(')(«) ~ • (3.18) 

n=0 n=0 


One may further determine the coefficients {a^,} and {bn} recursively by inserting these 
expansions into eq. (3.6) and collecting like powers of 7 and (1 —m). Without loss of 
generality, one may set oq = 1 and 5o = 0 . 

The values of frequency, for which the solution also satishes the Dirichlet condition 
at the boundary (for a fixed q), may also be expanded in powers of 7 , 

= + . (3.19) 


If the expansion (3.18) is truncated at some upper limit N and used throughout the 
computational domain 0 < m < 1, then the Dirichlet condition that 4/(0) vanish (for 
all 7 ) reduces to a set of algebraic equations 

N N 

«n(d;(°)) = 0, = 0, (3.20) 

n=0 n=0 
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Figure 1. The first few QNM frequencies, divided by 27rT, of the electromagnetic current 
operator for g = 0 (left) and q = 1 (right), evaluated at A = oo (red squares) and A = 
1000 (blue circles). The A = 1000 results include the 0 ( 7 ) corrections, but no higher order 
contributions. Lines have been inserted merely to guide the eye. 

whose roots yield (approximations to) and Carrying out this procedure for 
values of N sufficiently large that the results are stable turns out to be a viable com¬ 
putational strategy [51]. Our results for the hrst few QNMs are displayed in £g. 1 and 
reported in table 2 below, and fully agree with the findings of ref. [51]. 

We have also computed the photoemission spectrum (3.11) by solving the 7 - 
corrected equation for the transverse electric held, as in ref. [46], and then evaluated 
the hrst few moments ( 2 . 2 ) of the spectrum, obtaining the results shown in table 1 . 

As an alternative to the above approach, in which the QNM frequencies and mode 
functions are explicitly expanded in powers of 7 , one may directly solve the Schrodinger 
equation (3.6) with the potential (3.16) evaluated at some chosen value of 7 . Spectral 
methods provide an efficient numerical approach [58]. We write 

^(u) = (!-«)’'$(«), (3.21) 

so that the function $(«) is regular at both the horizon and boundary. In terms of <h, 
the explicit 7 -corrected QNM equation takes the form 

—u{l+u){l—u‘^) ^"{u) + u(1-|-m)^(1 — id}) ^'{u) + K{u) <h(u) = 0 , (3.22) 

with 
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K{u) = (1 - CP) + \u{l - 3 Cj‘^)) - \u^{l + CP) 


+ n4(^+“) 


11700 m + 2098482 _ 4752055 m® + 1838319 m^ 

+ f (4770 + 11700 M^ - 953781 m^ + 1011173 m®) 

+ (-4770 - 28170 + 1199223 m^) 1 . (3.23) 


Next, we expand 4> in a trnncated series of Chebyshev polynomials, 

N 

^u) = J2cnTn{2u-l), (3.24) 

n=0 

with Tn{z) = cos(m cos“^ ; 2 ).‘^ Reqniring eqnation (3.22) to be satisfied at the points 
Uj = I [1 — cos(j 7 r/A^)], j = 0,1,---N, which comprise a Chebyshev-Ganss-Lobatto 
grid, converts the differential eqnation (3.22) into a hnite matrix eqnation of the form 
A • M = 0. Here, m is a vector consisting of the iV+1 coefficients {c^} of the Chebyshev 
expansion, while A denotes a matrix whose entries on the j’th row are obtained by eval- 
nating eq. (3.22) at the j’th grid point.® Nonvanishing solntions of this homogeneons 
set of eqnations only exist when 

det(H)=0. (3.25) 

This determinant is a polynomial in a), whose roots rapidly converge (for hxed 

k) as the nnmber of grid points N increases. Evalnating these roots nnmerically is 
straightforward (for relatively modest valnes of N) given specihc valnes of the wavevec- 
tor q and the parameter 7 . 

Solving the QNM eqnation in the fashion described above yields valnes for the 
qnasinormal mode freqnencies with nonlinear dependence on 7 . One has, in effect, 
resnmmed a snbset of higher order contribntions to the QNM freqnencies which arise 
solely from the hrst order, i.e. correction to the snpergravity action. One may 

hope — althongh there is no gnarantee — that this is the dominant sonrce of all higher 
order contribntions. 

In hgnre 2 and table 2, we display the effects of the resnmmation for the hrst 
few QNM freqnencies, again at A = 1000 and both q = 0 and q = 1, and compare 
the resnlts to the previous unresummed values. As can be readily verihed, the size of 
the 0 ( 7 ) correction increases rapidly with the QNM mode number k, asymptotically 
growing like k‘^. This rehects the fact that hnite coupling corrections arise from higher 

^The Chebyshev polynomials {T„(z)} form an orthogonal basis in the Hilbert space with inner 
product {f,g) = dz f{z) g{z)/Vl - z^. 

®The row of this matrix corresponding to the m = 0 endpoint automatically enforces the Dirichlet 
boundary condition <i)(0) = 0. 
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Figure 2. The first few QNM frequencies, divided by 2ttT, of the electromagnetic current 
operator for g = 0 (left figure) or g = 1 (right figure) at A = 1000. Results obtained by 
directly solving the QNM equation at this value of A using spectral methods are shown as 
brown diamonds, while the red squares and blue circles show the same zeroth and first order 
results, respectively, previously displayed in fig. 1. Again, lines merely serve to guide the eye. 



u^iq = 0) 

or (? = 1) 

k 

0(7°) 

0(7^) 

resummed 

0(7°) 

0(7^) 

resummed 

1 

1-i 

1.073 - 1.005i 

1.068 -0.990i 

1.547- 0.8502 

1.558 - 0.8282 

1.557-0.8292 

2 

2-2i 

2.637- 1.797i 

2.237 - 1.794i 

2.399 - 1.874i 

2.525 - 1.6452 

2.477 - 1.722i 

3 

3-3i 

5.536 - 1.692f 

3.403 -2.55R 

3.323 - 2.8952 

3.957- 1.7912 

3.544 - 2.5182 

4 

A-Ai 

11.07+ 0.47i 

4.57 - 3.34z 

4.28 - 3.912 

6.37- 0.392 

4.67-3.312 


Table 2. The first four QNM frequencies, divided by 2ttT, of the electromagnetic current 
operator {w™}, /c = 1, • • • 4, for g = 0 (left) and q = 1 (right), at A = 1000. Respective 
columns show the results from the zeroth order, first order, and resummed approximations 
discussed in the text. 

dimension operators in the supergravity action (3.14). Although the size of the 0 ( 7 ) 
correction to the hrst QNM frequency appears modest at A = 1000, as seen in the 
top row of table 2, the 0 ( 7 ) correction exceeds the leading X = 00 value of the hrst 
QNM frequency at A = 139 (for g = 0), or A = 57.1 (for g = 1), above our benchmark 
phenomenological range. For the second mode (as noted previously) and all higher 
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-^breakdown (? — 0 ) 

•^breakdown 1 ) 

k 

0(7^) 

resummed 

0(7b 

resummed 

1 

139 

^ 22.5 

57.1 

< 5 

2 

382.2 

21.7 

195 

8 

3 

767.5 

21.3 

437 

14 

4 

1298 

^ 21.1 

793 

cs 16 


Table 3. Values of A below which the deviation of the QNM frequency w™ from its A = oo 
limit exceeds the A = oo value. Respective columns show the results obtained using either 
the first order or resummed approximations for the QNM frequency. 

modes, the “breakdown” values of A, below which the 0(j) correction exceeds the 
leading term, are even larger. 

In contrast, our resummed approximation for the QNM frequencies yields results 
which deviate from the A = oo values substantially less. As a concrete measure of this, 
table 3 compares, for the hrst few modes, the breakdown values of A below which the 
deviation of the QNM frequency from its A = cx) value exceeds the A = oo result, in 
either scheme.® For the resummed approximation, we hnd that these nominal break¬ 
down values of A are substantially smaller than the breakdown values of the hrst order 
results. For the hrst four modes, the breakdown values of the resummed approxima¬ 
tion lie within or below our benchmark phenomenological range of 10-40. Moreover, 
the breakdown values of the resummed approximation grow far less rapidly with mode 
number than do the 0(j) results. 

One may also apply our partial resummation scheme when evaluating the zero- 
frequency slope of the correlator which determines the electric conductivity of the 

= 4 SYM plasma via eq. (3.10). However, since the conductivity receives a much 
smaller 0 ( 7 ) correction than do the QNM frequencies, the ehect of the resummation 
is much more modest for the conductivity. At, for example, A = 1000 the relative size 
of the 0 ( 7 ) correction to the conductivity is about 8 x 10 “^, and our resummation 
increases this deviation from the A = 00 value by a further 5 x 10“®. In contrast to 
the situation with QNM frequencies, the nominal breakdown value of the resummed 
approximation for the conductivity is somewhat larger than the breakdown value of 
the hrst order approximation (A ~ 57.9 instead of 39.7). 

®When computing resummed approximations, the convergence of both the Frobenius (3.18) and 
spectral (3.24) expansions progressively degrade as the value of 7 is increased, making precise deter¬ 
minations of breakdown values of A challenging. 
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3.2 Stress-energy correlator 


The dynamics of linearized metric pertnrbations determine the stress-energy tensor 
correlator. We will focns on the £ = 2 or shear channel, for which it is snfficient to 
consider as the only non-zero component of the perturbation when the wavevector 
points in the z-direction. ' After Fourier transforming with respect to the boundary 
coordinates, the rescaled perturbation 

11 

Z{u) = ^ 5g^y{u) (3.26) 

satishes the 0 ( 7 ) corrected equation of motion [59], 

- Z''{u) + P{u) Z\u) + Q{u) Z{u) = 0 , (3.27) 


with 


and 


1 -I- 


P(u) = + I 7 (600 u - 2306 - 3171 - 3840 f u^) (3.28) 

ut[u) 

Q{u) = -(I + 3 O 7 ) 


- q^f{u) 


7 


4/(«) 


uf{uy 

50m2- 275m® + 225m^ 

+ (600 u - 2856 + 2136 u^) 

+ f (-300 u + 3456 - 6560 + 3404 u’^) 

+ (768m^ + 768m®) 


(3.29) 


Boundary conditions for hnding quasinormal modes are the same as discussed ear¬ 
lier: infalling behavior at the horizon and Dirichlet at the boundary. Solutions up to 
0 ( 7 ) of the above equation were found in ref. [50] using the Frobenius expansion tech¬ 
nique discussed in the previous subsection. One may, however, also solve the equation 
directly for specihc values of 7 , just as we did for the electromagnetic current corre¬ 
lator, using spectral methods. Figure 3 shows a comparison of the results of the two 
approaches for A = 500 and g = 0, with table 4 listing explicit values. Once again, 
at this value of A (which is still well above the phenomenologically relevant range) we 

^Note that in ref. [57], the channel with £ = 2 rotational symmetry about the wavevector was 
referred to as the “scalar channel” because the corresponding metric perturbation satisfies the same 
equation as a minimally coupled massless scalar, and the t = 1 or vector channel was referred to as the 
“shear channel”. Our terminology is motivated by the fact that the shear viscosity can be obtained 
from the zero-frequency limit of the correlator in the £ = 2 channel. 
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fHe(w) 

Figure 3. The first few QNM frequencies, divided by 27rT, for the shear channel of the 
stress-energy correlator, evaluated for g = 0 and A = 500. Results obtained by directly 
solving the QNM equation at this value of A using spectral methods are shown as brown 
diamonds, while the red squares and blue circles show results truncated at zeroth and hrst 
order in 7 , respectively. As before, lines merely serve to guide the eye. 



^shearj-^ = Q) 

k 

0(7°) 

o(y) 

resummed 

1 

1.560 - 1.373i 

1.581 - 1.3562 

1.579 - 1.3562 

2 

2.585 - 2.382i 

2.723 - 2.2532 

2.673 - 2.2772 

3 

3.594 - 3.385i 

4.093 - 2.8992 

3.789 - 3.1332 

4 

4.60 - 4.392 

5.92 - 3.072 

4.91 - 3.982 


Table 4. The first four QNM frequencies, divided by 2'kT, fe = 1, • • • 4, in the shear 

channel of the stress-energy correlator, for q = 0 and A = 500. Respective columns show 
the results from the zeroth order, first order, and resummed approximations discussed in the 
text. 

observe a substantial difference between the 0(7) results and our resummed values, 
with the resummation decreasing the difference from the A = cx) limit. In parallel with 
the previous subsection, we show in table 5 the breakdown values of A, computed with 
both hrst order and resummed approximations, below which the hrst few shear channel 
QNM frequencies deviate from their A = 00 limits by more than their X = 00 values. 
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-^breakdown 

k 

0(7b 

resummed 

1 

27.6 

< 1 

2 

71.2 

< 1 

3 

135.5 

< 1 

4 

221 

< 1 


Table 5. Values of A below which the deviation of the QNM frequency from its A = oo 

limit exceeds the A = oo value. Respective columns show the results obtained using either 
the hrst order or resummed approximations for the QNM frequency. 


We find that for the shear channel qnasinormal freqnencies, onr resnmmation scheme 
leads to nominal breakdown valnes of A for all modes np to fc = 4 lying below our 
benchmark phenomenological range of 10-40. 

Similarly to the plasma conductivity considered earlier, the shear viscosity is de¬ 
termined by the zero-frequency slope of the retarded correlator of T^y evaluated at 
vanishing wavenumber, 

h = - lim Im- {u, 0). (3.30) 

Including 0 ( 7 ) corrections, it can be shown that this correlator is given by [43] 


G 


ret 

xy^xy 




{rir Z\u) 

u^o 47r2 uZ{u)' 


(3.31) 


where Z{u) is a solution to eqs. (3.27)-(3.29) at g = 0 that satishes infalling boundary 
conditions at the horizon, and = vrTL ^/(1 -|-15 7 ) is the horizon position of the 7 = 0 
geometry.® Using the Frobenius method to solve for the metric perturbation to linear 
order in 7 now leads to 


G 


ret 

xy,xy 


{oj, 0 ) = —i 


Stt^T 


(1 + 195 7) + C>(u;^72), 


(3.32) 


from which one can extract the result found in ref. [44], 7 = | {1 + 135 7 -|- • • •). 

One may also apply our partial resnmmation scheme when evaluating the zero- 
frequency slope of the correlator (3.31). As with the conductivity, since the shear 
viscosity receives a much smaller 0 ( 7 ) correction than do the QNM frequencies, the 

®This 7 = 0 horizon position, used consistently in ref. [43], differs from the horizon position = 
'kTL'^/{I -I- ^ 7 ) of the 7 -corrected metric derived in ref. [60] and used throughout refs. [19, 45, 46]. 
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effect of the resummation is much more modest for the this transport coefficient. At, 
for example, A = 1000 the relative size of the 0(j) correction to the viscosity is about 
6 X 10“"^, and our resummation increases this deviation from the A = cx) value by a 
mere 10“^. As with conductivity, the resummed approximation for the shear viscosity 
leads to a somewhat larger nominal breakdown value of A as compared to the hrst order 
approximation (A ~ 14 instead of ~ 7). 

4 Conclusions 

Our examination of £nite-A corrections to holographic results for thermal quantities 
is motivated by an obvious desire to understand more clearly the applicability of 
gauge/gravity duality to the physics of quark-gluon plasma as produced in real heavy 
ion collisions. In particular, how large an error is made when the system is mod¬ 
eled as inhnitely strongly coupled, as is customary in most holographic calculations of 
non-equilibrium dynamics performed in the supergravity limit? We approached this 
problem by comparing the behaviors of the strong coupling expansions for a variety of 
thermal quantities, for which the leading order finite coupling correction, of order A“^/^ 
in the’t Hooft coupling, is known. Our comparison has singled out quasinormal mode 
frequencies as quantities for which the hnite coupling corrections appear particularly 
problematic at phenomenologically interesting values of the’t Hooft coupling. We dis¬ 
covered, however, that a partial inclusion of higher order contributions, generated by 
the leading corrections to the supergravity action, leads to a dramatic reduction in the 
predicted size of hnite coupling effects in quasinormal mode frequencies. 

One clear, but unsurprising, message of our comparison is that notions of strong 
(or weak) coupling domains can depend rather strongly on the physics observable of 
interest. For some quantities, ’t Hooft couplings of order 10-40 may be easily accessible 
via a one or two term strong coupling expansion, while for other quantities this is clearly 
not the case. The full reasons underlying such behavior are unclear, but some insight 
may perhaps be drawn from similar issues at weak coupling. There, an extensive 
amount of work has been devoted to the calculation of high order perturbative results 
for a number of equilibrium thermodynamic quantities. An important lesson from 
this work is that the convergence of perturbation theory is intimately related to the 
relative inhuence of different momentum or energy scales contributing to the observable 
in question. Weak coupling expansions for physical quantities that are dominantly 
sensitive to the hard thermal scale of 2ttT are found to behave signihcantly better than 
expansions of quantities which are more sensitive to the soft gT and ultrasoft g^T scales 
which originate from electro- and magnetostatic screening, respectively. The degraded 
perturbative stability arises from the dependence on the ratio of these scales, which 
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appears in the form of contributions suppressed only by single powers (and logarithms) 
of g, instead of as = g‘^/{4:7i). Ultimately, this reflects the diverging infrared sensitivity 
of the Bose-Einstein distribution of gluonic helds. 

Whether it is possible to interpret differences in stability of strong coupling ex¬ 
pansions in an analogous manner remains to be seen. However, considering the fact 
that the 0{a'^) corrections to A = cxd results originate from higher derivative operators 
added to the supergravity action, and that the leading hnite coupling corrections to the 
QNM frequencies can be seen to grow rapidly with the mode number, perhaps such ex¬ 
pectations are not completely unfounded. In the meantime, trying to generalize many 
more holographic calculations to include at least the leading hnite ’t Hooft couplings 
would clearly be worthwhile. 
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